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Abstract. Large transverse momentum distributions of identified particles observed at RHIC are analyzed 
by a relativistic stochastic model in the three dimensional (non-Euclidean) rapidity space. A distribution 
function obtained from the model is Gaussian-like in radial rapidity. It can well describe observed transverse 
momentum pr distributions. Estimation of radial flow is made from the analysis of pr distributions for 
p in Au + Au Collisions. Temperatures are estimated from observed large px distributions under the 
assumption that the distribution function approaches to the Maxwell-Boltzmann distribution in the lower 
momentum limit. Power-law behavior of large pr distribution is also derived from the model. 



PACS. 25.75.-q Relativistic heavy ion collisions 
20.50.Ey Stochastic models 



13.85.Ni Inclusive production with identified hadrons 



1 Introduction 

At RHIC colliding energy of nuclei becomes up to 200 
AGeV, and thousands of particles are produced per event. 
To describe such many particle system, a sort of collective 
approach will be useful. One-particle rapidity or pseudo- 
rapidity distributions observed at RHIC are well described 
by the Ornstein-Uhlenbeck process 0121 • 

In order to treat transverse momentum px distribu- 
tions, we can extend one dimensional stochastic equation 
to a three-dimensional one using variables, longitudinal ra- 
pidity and two transverse momenta. The fundamental so- 
lution is Gaussian in these variables. Even if the formula 
is integrated over the azimuthal angle in the transverse 
momentum plane, the distribution is only slightly modi- 
fied from the Gaussian in the low px region. However, the 
observed pr distributions at RHIC have long tails com- 
pared with an exponential distribution, and cannot be de- 
scribed by the formula. Therefore relativistic treatment 
for the transverse direction would be needed as long as 
the stochastic approach is adopted. 

In reference [H], an empirical formula for large px dis- 
tributions at polar angle 9 — 7r/2, 
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was proposed from the analogy of Landau's hydrodynam- 
ical model. Polar angle 9 (0 < 9 < tt) is measured from 
the beam direction of colliding nuclei or incident particles. 
In equation (JJJ, E denotes energy of an observed particle, 
Lt is a parameter, and y-r is called the "transverse rapid- 
ity" . Equation well describes large px distributions for 
p + p — > 7T° + X and p + p — ► 7r ± + X. Recently, large 
Pt distributions for tt° in pp collisions observed at RHIC 
are also analyzed by QJ 4 . However, it cannot be derived 
from the hydrodynamical model. 

The transverse rapidity is defined in the geodesic cylin- 
drical coordinate system in the three dimensional rapidity 
space, where longitudinal rapidity y, transverse rapidity £ 
and azimuthal angle are used. The longitudinal rapidity 
y, and the transverse rapidity £ are defined respectively 
as, 



i E+ P L c 
rriT 



In 



IPt] 



where E,pl, Pt and m denote energy, longitudinal mo- 
mentum, transverse momentum, and mass of the observed 
particle, respectively, and m-r = */ pE + to 2 . It should be 
noted that yx coincides with £, only if 9 = tt/2, namely 

PL = 0. 

As for the relativistic approach to stochastic equation, 
we consider the diffusion equation in the three dimen- 
sional rapidity space or Lobachevsky space, which is non- 
Euclidean. In order to classify the longitudinal and trans- 
verse expansion, it would be appropriate to consider the 
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diffusion equation in the geodesic cylindrical coordinate 

D 



df D 8 2 f D d ( . , . , , 

"tt - = o — rr^r H — . : t r— -7— sinn t cosh E 

dt cosh 2 £ cV sinh £ cosh £ <9£ V 

D d 2 f 
+ sinh 2 ^0 2 ' 



(2) 



where D denotes a diffusion constant. However, we cannot 
solve at present, as far as we are aware. Therefore, we 
should consider somewhat simpler case. 

We have proposed the relativistic diffusion model where 
an effect of radial flow is not included, and analyzed large 
pr distributions for charged particles [H] and identified 
particles 6 in Au+Au collisions. The distribution function 
of it is Gaussian-like in radial rapidity p, which coincides 
with £ at 8 — tt/2, and resembles with (JJJ at 9 — ir/2. 

In section 2, the relativistic diffusion model, where the 
flow effect is taken into account, is briefly explained. In 
section 3, large pr distributions in p + p — > tt° + X and 
in Au + Au — > 7T° + X at ^/snn — 200 GeV observed at 
RHIC 7,8,9 are analyzed. In section 4, a flow effect or a 
transverse flow velocity is estimated from px distributions 
in Au + Au — ► p + X , and pt distributions for 7r°, 7r~ and 
X~ are analyzed by the use of the transverse flow velocity 
estimated from px distributions of p. A possible interpre- 
tation of the dispersion of the Gaussian-like distribution 
is shown in section 5. The higher transverse momentum 
limit of our model is discussed in section 6. Final section 
is devoted to summary and discussions. 



2 Diffusion equation in the three dimensional 
rapidity space 

For simplicity, we consider the diffusion equation with ra- 
dial symmetry in the three dimensional geodesic polar co- 
ordinate system, 



df d a 



2 2i 

"dp 



with an initial condition. 



f(p,t = 0) 
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In equation (J3J), p denotes the radial rapidity, which is 
written with energy E, momentum p and mass m of ob- 
served particle, 



p = In 



E+\p\ 



(5) 



Inversely, energy and momentum are written respectively 
as 



E = mcoshp, |p| = \lp\ + p T = msmhp. (6) 



The solution rTU] of © with the initial condition (TJJ 
is given by 



f(p,p ,t) = 
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A physical picture described by 10 and (TJJ is as fol- 
lows; After a collision of nuclei, particles are produced 
at radial rapidity po expressed by (TJJ). Then those par- 
ticles propagate according to the diffusion equation @. 
In the course of the space time development, energy is 
supplied from the leading particle system (collided nu- 
clei) to the produced particle system. Number density of 
particles becomes lower and at some (critical) density, in- 
teractions among secondary particles cease, and particles 
become free. 

From equation (JJJ, the following equation is ob- 
tained; 



f{p,t) = lim . f(p,p ,t) 

po— >0 



= (47rDi) 
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Transverse momentum (rapidity) distributions at fixed po- 
lar angle 9 are analyzed by 10 or iJSJ, where transverse 
momentum is given by |pr| = m sinh p sin 9. 

In order to clarify our approach to the previous non- 
relativistic approach 12,13 , equation Q is compared with 
the solution for the radial symmetric diffusion equation in 
the Euclidean space in appendix A. 



3 Analysis of p T distributions observed at 
RHIC 

At first, transverse momentum (rapidity) distributions of 
identified particles observed by the PHENIX collabora- 
tion 0IH1OI1 are analyzed by the use of (JSJ) , 

f(p,t)=C(2na(t) 2 y 3/2 e-^ 2 / 2 
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with parameters, C and <x(t) 2 1 . 

As the data are taken at 8 = it/2 or p^ — 0, the 
identity p = £ = yr is satisfied. In the analysis, we use 
m^o = 135, m^- = 140, m K - = 494, and m p = 938 MeV. 

The result on the px distribution in p + p — > ir° + X is 
shown in Fig.f^ Solid curve in Fig.fjis drawn by the use 
of (PJ, parameters of which are estimated with the least 
square method, and are shown in Table ^ Observed pt 



1 Somewhat different analyses based on statistical models 
and the different stochastic model have been made in |14) 
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Fig. 1. pt distribution for p + p — > 7r° + X at y — 7 
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Table 1. Parameters on px distributions estimated by © in 
p + p -> 7T° + X at y = at V» = 200 GeV 



<x(t) 2 



X 2 /n.d.f. 



13547.3 ± 941.7 0.603± 0.004 12.3/15 



Table 2. Parameters on pt distributions estimated by in 
+ Au -> 7T° + X at j/ = at ys^zv = 200 GeV [HI 



centrality 


C 


oitf 


xVn-d.f 


min.bias 


45066.0 ± 4542.7 


0.586 ± 0.004 


28.0/22 


00-10% 


146819.0 ± 16095.3 


0.574 ± 0.005 


27.3/16 


10-20% 


103772.0 ± 11348.3 


0.580 ± 0.005 


20.0/15 


20-30% 


69568.0 ± 7647.8 


0.586 ± 0.005 


14.9/13 


30-40% 


43168.0 ± 4870.7 


0.592 ± 0.006 


15.7/13 


40-50% 


25292.0 ± 2955.2 


0.600 ± 0.006 


9.4/12 


50-60% 


13800.3 ± 1629.2 


0.600 ± 0.007 


12.4/12 


60-70% 


5477.2 ± 699.8 


0.617 ± 0.008 


8.2/10 


70-80% 


2498.80 ± 336.74 


0.616 ± 0.009 


7.2/9 


80-92% 


1275.89 ± 216.03 


0.610 ± 0.011 


5.0/8 



distributions for ir° in p + p collisions are well described 
by © from 1.2 GeV/c to 13.2 GeV/c. 

The results on pt distributions in Au + Au — * tt° + X 
are shown in Fig.[3]and Table|3 As can be seen from Fig. [2] 
and Table |3 observed p^ distributions on ir° in Au + Au 
collisions are well described by I© from centrality 0-10% 
to centrality 80-92 %. 



Fig. 2. p T distribution for Au + Au -*• 7r° + X at y = [H] 



4 Estimation of radial flow 

We have analyzed large pt distributions in Au + Au — > 
n°+X in the previous section. However, if the effect of flow 
in the transverse direction is taken into account, estimated 
values of parameters would be different from those shown 
in Table □ 

The effect of transverse flow in Au + Au collisions at 
\J s nn — 200 GeV is estimated from the observed pr dis- 
tribution in Au + Au — > p + X from pr = 0.65 GeV/c to 
Pt = 4.25 GeV/c by the use of the formula, 



f(p,Po,t) = 



C 



2iry/2ir<r(t) 2 
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sinh 



sinh po sinh p 



exp 



pl 
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(10) 



The results are shown in Fig.|21and Table|3J From TablcEl 
we can obtain the transverse flow velocity vq using the 
relation vo — tanhpo- It varies from vq = tanh0.82 = 
0.68 at centrality 00-05% to v = 0.42 at centrality 80- 
92% (v = 0.39 at centrality 70-80%). As the relation of 
centrality to the average number of participant nucleons 
(Nprt) is given in Refs. |Zj and [H], we can fit the estimated 
values of po as a function of (N prt ) ; 



po = 0.113 ln(a:) + 0.164, x = (N prt ). 



(11) 



The result is shown in Fig. fj] 

Then, we analyze observed pt distributions with po 
given by (jllfl . and (N prt ), which is shown in The re- 
sults on 7T°, tt~ and K~ are shown respectively in TablcsQl 
0and El 
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Table 3. Parameters on pr distributions estimated by ill (It in Au + Au -^p + Xaty = 0at ^/sjvjv = 200 GeV [5] 



centrality C °~(i) 2 Po X 2 / n -d-f 



00- 


05% 


17.713 ± 0.626 





138± 





008 





820± 





029 


4 


5/19 


05- 


10% 


14.974 ± 0.538 





141± 





009 





803± 





030 


5 


4/19 


10- 


15% 


12.508 ± 0.454 





144± 





009 





795± 





032 


4 


4/19 


15- 


20% 


10.492 ± 0.389 





147± 





010 





778± 





034 


3 


7/19 


20- 


30% 


7.762 ± 0.293 





148± 





009 





759± 





034 


3 


2/19 


30- 


40% 


5.175 ± 0.209 





158± 





011 





705± 





041 


1 


0/10 


40- 


50% 


3.168 ± 0.137 





167± 





001 





647± 





051 


2 


3/19 


50- 


60% 


1.8169 ± 0.0877 





172± 





001 





587± 





066 


1 


8/19 


60- 


70% 


0.8949 ± 0.0489 





191± 





001 





459± 





127 


1 


0/19 


70- 


80% 


0.3638 ± 0.0239 





190± 





001 





409± 





211 


6 


1/19 


80- 


92% 


0.1821 ± 0.0157 





167± 





001 





444± 





195 


2 


6/17 



Table 4. Parameters on pr distributions estimated by (11 01 with il l 1 1 in Au + Au — ► ir° + X at y = at -^/sjvjv = 200 GeV 



centrality 




C 




a 


(tf 


Po 


X 2 /n.d.f 


min.bias 


27977.0 


± 


2764.2 


0.500 


± 


0.004 


0.694 


29.8/22 


00-10% 


72773.0 


± 


7625.2 


0.469 


± 


0.005 


0.818 


24.2/16 


10-20% 


56883.0 


± 


5936.0 


0.479 


± 


0.005 


0.781 


17.0/15 


20-30% 


41647.0 


± 


4386.1 


0.490 


± 


0.005 


0.742 


10.3/13 


30-40% 


27513.0 


± 


3012.5 


0.503 


± 


0.005 


0.699 


13.1/13 


40-50% 


17515.0 


± 


2002.8 


0.518 


± 


0.006 


0.651 


7.6/12 


50-60% 


10559.0 


± 


1224.0 


0.526 


± 


0.006 


0.595 


9.2/12 


60-70% 


4548.40 


± 


578.86 


0.553 


± 


0.008 


0.531 


7.2/10 


70-80% 


2241.30 


± 


313.24 


0.564 


± 


0.009 


0.457 


6.4/9 


80-92% 


1178.50 


± 


204.72 


0.574 


± 


0.011 


0.372 


4.9/8 



Table 5. Parameters on pr distributions estimated by (11 01 wit 1)11^ in Au + Au — > n + X at y — at v / sjvjv = 200 GeV 



centrality C °"(^) 2 Po x 2 / n -d-f 



00- 


05% 


32150.0 ± 846.3 





533± 





003 





826 


357.4/26 


05- 


10% 


27161.0 ± 718.2 





546± 





003 





808 


308.0/26 


10- 


15% 


22407.0 ± 594.1 





559± 





004 





790 


271.1/26 


15- 


20% 


18795.0 ± 501.4 





566± 





004 





781 


217.5/26 


20- 


30% 


14332.0 ± 383.4 





582± 





004 





742 


195.3/26 


30- 


40% 


9337.0 ± 254.4 





603± 





004 





699 


136.9/26 


40- 


50% 


5773.6 ± 161.0 





617± 





004 





651 


108.6/26 


50- 


60% 


3365.10 ± 96.36 





629± 





005 





595 


91.0/26 


60- 


70% 


1729.70 ± 51.22 





640± 





005 





531 


80.8/26 


70- 


80% 


784.01 ± 24.37 





655± 





006 





457 


54.4/26 


80- 


92% 


393.59 ± 12.88 





664± 





007 





372 


39.9/26 



As can be seen from Tables 0] and [5l fitting of our 
calculation by (|ll)f) to the data is better than that by 10 
except for the minimum bias events in Au + Au — > tt° + X. 
The fit to 7r~ distributions is not so good as that to ir° 
distributions. The result at lower centrality is no good, 
but it is much improved as the centrality increases. 



The results by UlUjl on it and K are also better than 
those by J^J, which are not shown for brevity. 



5 Possible interpretation of parameter a(t)' 



In order to get a possible interpretation of parameter a(t) 2 , 
an approximate expression of @ in the small \pr \ region 
is taken. When p << 1, |p| = msinhp ~ mp. Then, © 
reduces to 



f(p, t) ~ exp 



2a(tX< 



(12) 
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Table 6. Parameters on pr distributions estimated by il l with il l 1 1 in Au + Au — » K~ ± X at y = at ^/sjvjv = 200 GeV [5] 



centrality 


C 


a(t) 


pa 


X /n.d.t 


00-05% 


254.63 ± 8.93 


0.258± 0.008 


0.826 


0.52/14 


05-10% 


212.58 ± 7.45 


0.266± 0.008 


0.808 


0.62/14 


10-15% 


175.59 ± 6.20 


0.269± 0.008 


0.790 


0.87/14 


15-20% 


145.12 ± 5.17 


0.271± 0.008 


0.781 


0.92/14 


20-30% 


108.69 ± 3.89 


0.286± 0.008 


0.742 


1.07/14 


30-40% 


70.370 ±2.579 


0.290± 0.008 


0.699 


1.70/14 


40-50% 


42.070 ±1.577 


0.302± 0.008 


0.651 


1.80/14 


50-60% 


23.404 ±0.918 


0.311± 0.009 


0.595 


2.76/14 


60-70% 


11.308 ±0.465 


0.320± 0.009 


0.531 


3.51/14 


70-80% 


4.782 ±0.215 


0.329± 0.011 


0.457 


3.75/14 


80-92% 


2.274 ±0.113 


0.342± 0.013 


0.372 


7.59/14 
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Fig. 3. pr distribution for Au + Au 
by (IIOII . where po is included. 



• p at y = 0, analyzed 



Fig. 4. Dependence of po to average number of participants 
in Au ± Au collisions estimated by Hill 



If we assume that equation i|12|) should coincide with the 
Maxwell-Boltzmann distribution, 



/(v) = 



\2TTk B T y 



3/2 




mv 2 " 


J exp 


. 2k B Tl 



where k B is the Boltzmann constant, then we have an 
identity, 



k B T = ma(ty 



(13) 



Equation (|13f) is obtained in the lower momentum limit, 
namely, with the condition that p « 1. 

From l|13|) . we can estimate the temperature k B T of 
inclusive reactions for observed particle with mass m. The 
results are shown in Fig. The estimated temperature 



k B T from the pr distributions for 7r° is from 63.3 MeV at 
centrality 5-10% to 77.5 MeV at centrality from 80-92%, 
and that for tt~ is from 74.8 MeV at centrality 0-5% to 
93.1 MeV at centrality from 80-92%. The temperature of 
7T° or 7r~, estimated from H13(l . is somewhat lower than 
that from the exponential function of px- This would be 
affected by the assumption p << 1 used to obtain the 
identity (fT3^l . 

The estimated temperature from the px distributions 
for K- is from 128.0 MeV at centrality 0-5% to 163.5 
MeV at centrality from 80-92%. That for p distributions 
is from 128.5 MeV at centrality 0-5% to 172.6 MeV at cen- 
trality from 80-92%. The results on jj distributions shown 
in Fig. are recalculated by the use of ijT2|l . 
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for any positive number e. Therefore, we can approximate 
p in the exponent in i|14|) as constant, which is written by 
2co, within some finite transverse momentum range. Then 
(|10[l is reduced to the form, 



f(p,Po,t) 



rn 



(eo-po)Mt) 2 + l 



m T + Px , 

(oo-po)/ff(t) 2 +l 



(15) 



From equation (|15|l , one can see that equation (|10|l shows 
power-law behavior in the higher transverse momentum 
limit, and that the power becomes smaller as po increases, 
and it also becomes smaller as u(t) 2 , which should increase 
with the colliding energy ^/sjvjv, increases. 



Fig. 5. Centrality dependence of temperature ksT estimated 
by HUH with 1111 from px distributions in Au + Au collisions 

6 Higher transverse momentum limit of the 
model 

In the 1970's, when large px distributions are observed in 
accelerator experiments, many models, which have power- 
law behavior in pt are proposed |15| . In |16j , a model for 
Pt distribution, inspired by QCD, is proposed as {po/(px+ 
Po)) n . In the following, |px| is abbreviated to px- It ap- 
proaches an exponential distribution of pt for px — > 0, 
and a power function of px for px — * oo. 

In this section, it is derived that equation Q at 8 = 
tt/2 shows power-law behavior in px in the higher px limit, 
when the identity p = 1ii((tot + px)/m) is used. 

The following two terms in (J7J can be rewritten in the 
power form of tot + Pt as; 



sinh 



PoP 
<7*(i) 



1 _ e -2p p/a 2 (t) 



.Pop/ff it) 



poM*) 



^ 1 / m T + px 
~ 2 V to 
1 - e~ 2p If m T + px 



sinh p = 



Radial rapidity contained in the Gaussian-like part in I|1U|) 
is rewritten as, 



exp 



P 



2a(t) 2 



m T + Px 



rn 



-p/2a(t) 2 



(14) 



As is well known that a logarithmic distribution of x has 
weaker dependence compared with a power-law distribu- 
tion of a; with any positive exponent in the limit of x — > oo, 
we have, 



lim — 

PT — >00 p T 



,. 1 , f mx+px 
lim — In 







7 Summary and discussions 



In order to analyze large px distributions of charged par- 
ticles observed at RHIC, the relativistic stochastic pro- 
cess in the three dimensional rapidity space, which is non- 
Euclidean, is introduced. The solution is Gaussian-like in 
radial rapidity, where the radial flow rapidity po is in- 
cluded. It is very similar to the formula proposed in |3] at 
9 = tt/2. 

Transverse momentum distributions for ir° in pp colli- 
sions at y = at y/s = 200 GeV observed by the PHENIX 
collaboration are analyzed. Observed px distributions are 
well described by the formula © from 1.2 GeV/c to 13.2 
GeV/c 

In Au + Au collisions, firstly, we analyzed the data 
for 7T° using ©, where flow effect is not included, i.e. 
po = 0. The formula © well describes the observed px 
distributions. 

Next, we have analyzed the observed px distributions 
for Au + Au —fp + X using (|10fl ■ In the formula, non-zero 
radial rapidity p at the initial stage, which corresponds 
to the radial flow velocity v = tanhp 0j is included. 

The temperature is estimated under the assumption 
that (0} approaches to the Maxwell-Boltzmann distribu- 
tion when p « 1 or px << to. 

In our analysis, non-zero value of the radial flow ra- 
pidity (or radial flow velocity) has an influence on the 
estimation of temperature. One can see the effect of po 
from (|15fl . In addition, all of the minimum chi-squared 
values by the use of i|10[l with l|ll|l become smaller than 
those with po = except for the minimum bias events in 

Au + Au — > 7T° + X. 

We have derived that our formula (|10fl shows power- 
law behavior in px in the higher transverse momentum 
limit. Therefore, it behaves like a Gaussian distribution in 
Px when px << to, and like a power-law distribution in 
Px when px >> 1 • 
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A Relation to non-relativistic approach 

The radial symmetric diffusion equation in the n (n = 2, 3) 
dimensional Euclidean space is written as, 



Of =D d 2 f | D(n-l) df 
dt dr 2 r dr 

The solution of Ijl6(l with the initial condition 

5(r - r ) 



/M = 0) 



2Tr(2r o y 



is given by [*T7). 



x (r Q r) 



4Dt 



-(n-2)/2 



z (n-2)/2 



/ r r \ 
\2Dt) 



(16) 



(17) 



(18) 



In the case of two dimensional expansion, n = 2, l|18fl 
coincides with the formula in 12 , if 2Dt — m/T, r = 
|Pt|/to, and r = |pto|/^- 
Using the identity 
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sinhz = y/TTz/2I 1 / 2 (z), 
we can rewrite {7J| as 



f{P, Po,*) 



1 



(POP) 



1/2 



8ttL>< 
x exp 



smb. po sinh p 



p 2 + pl 



ADt 



In the case of the three dimensional expansion in the Eu- 
clidean space, comparing l|19l) with (|18fl at n = 3, one can 
recognize that <|19|) is a relativistic extension of (|18f) . and 
that (|19|l reduces to Ijl8(l in the low momentum region, 
|p| << 1, where |p| = msinhp ~ mp and sinhpo — Po- 
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